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Abstract 

Recent exact analytical results developed for the random number generators with 
taps are reported. These results are applicable to a wide class of algorithms, includ- 
ing random walks, cluster algorithms, Ising models. Practical considerations on the 
improvement of the quality of random numbers are discussed as well. 
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1 Introduction 



The large-scale simulations in statistical physics use some deterministic pro- 
cedures to generate a sequence of uniformly distributed (pseudo)random num- 
bers. It is possible to generate 10^ numbers per second and 10^^ numbers in 
100 days on the best processors. The widely known linear congruential meth- 
ods cannot be used in a such simulations on 32-bit computers because of the 
short period T < 2^"^ ^ 10^ of the generated sequence (|^). The next most used 
class of generators are based on the linear recursion 



Xn = {Xn-g OP Xn-p) mod 771, 



where m is the word length and p and q are feedback taps p > q. Here OP 
stands for the exlusive-or operation for shift registers, + operation for the 
lagged Fibonacci (LF) and — operation for the substract-with-carry (SWC) 
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recipes. In the latter case one has to add unity if the result of the previous 
substraction was negative. The correlations in 2-tap generators were found to 
be crucial in large-scale simulations of Ising model, self-avoiding random walks 
(H) , percolation (^) and cluster algorithms (H) (see, also, (||)). It was shown 
recently by Heringa, Blote and author that these deviations can be rigorously 
estimated numerically and even analytically(|^). Using this idea, Blote and 
author discovered the scaling of systematic deviations for Ising model (|^. 

The method gives us some rigorous and quantitative way to estimate the 
possible level of deviations in the results of simulations. So, we propose some 
practical definition of goodness of a concrete recipe for generating random 
numbers. This definition is based on the concept of observability: the recipe 
should be considered as a good one if it passes all known tests using the 
sequence of the pseudo-random numbers generated in a couple of month on 
the nowadays computers. So, for today we need those recipes which do not 
produce observable deviations when 10^^ random numbers are used. 



2 Random walk test 



The test is based on the one-dimensional directed random walk (P): a walker 
starts at some site of an one-dimensional lattice and, at discrete times i, either 
he takes a step in a fixed direction with a probability /ij or he stops with a 
probability z/j = 1 — yUj. In the case of unbiased probabilities /ij = /i and the 
probability of a walk with length n is then P{n) = /i""^ (1 — /x). We could 
apply some particular recipe of random number generation to calculation of the 
probabilities /ij. Clearly, they would differ from because of the correlations 
in the recipe (1). Therefore, it was proposed in (H; ^ to calculate the relative 
deviations 

SP^n) - ^ - 1 (2) 



of the probability of the walk P*{n) = v* nr=i(/^i)* ^ from the unbiased prob- 
ability P{n). 

Some results for SR anf LF generators from and for SWC one from are 
presented in the Table 1. The first interesting fact is that the lagged Fibonacci 
and substract-with-carry leads to the same deviations up to the order of 2~™, 
which is usually of the order of 10~^. Next, the deviations for shift registers are 
of the same nature and and of the same order as for LF and SWC. Both these 
statements contradict to the widespread belief that LF and SWC produce 
"better" randomness. 
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Table 1 

Deviations 5P{n) of the probability of walk length n due to the correlations in 
the random number sequences. Corrections of the order of 2~™ are not included in 
results for LF and SWC. 



length n 


2-tap SR 


2-tap LF = SWC 


P 


M 


1-2^ 
2fi 


p+l 




(3m-1)2 . 
4^-1 ^ 



3 Scaling of the deviations 



We performed Wolff simulations of the 2D Ising model at criticality, using SR 
with feed-back positions (p, g) = (36,ll), (89,38), (127,64) and (250,103) for a 
variety of the lattice sizes L (0). We found that the data could be collapsed on 
a single curve with the appropriate rescaling of deviations, say, of the energy 
SE = p'^-^^SE and the system size L = p~^-^^(^)L. This enables us to bound 
deviations of energy E 



5E< 0.3 L-°-84 



(3) 



of the specific heat C 



-0.21 ^-0.42 



6C< 0.85 L-^ '^ p 



(4) 



and of the universal ratio Q = (m^)^/ (m^), where m is the magnetisation, 

5g< 0.244 L-°-^5p-°-^^ (5) 

The analysis has been performed also for 3D Ising model, and it yields 

~6C ^ (6) 



and 



5Q ^ /-o^ 



(7) 



These estimations are very useful in the planning of computer experiments 
([TT|). Their knowledge is equivalent to the knowledge of the accuracy of the 
experimental equipment. 
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4 Consclusion 



Using the method developed in (H) we found some support for the previ- 
ously proposed modifications of the 2-tapped methods. Correlations would be 
sharply diminished by one of the following tricks: 

i) concept of luxuries |^), i.e., generate r > p random numbers but use 
only p of them; 

ii) by decimation of the random number sequence, or, by using a combination 
of the two or more generators ([T0| ; |lll) ; 

iii) by using four-tap generators ( p^ ||). 



Our results are shown to be important for the simulations of Ising model, the 
percolation problem, for random walks and cluster algorithms. 

So, instead of the qualifications "bad" or "good" we propose to apply to 
random number generators a more detailed description of their features (^, 
for example, that no deviations would be observed in the simulation of the 
above mentioned problems when up to 10^^ random numbers are used. 

Discussions and collaboration with H. Blote, J. Herringa and P. Butera are 
kindly acknowledged. 
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